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ABSTRACT 

The need for unified constitutive models was recognized more than a decade ago in 
the results of phenomenological tests on monolithic metals that exhibited strong 
creep— plasticity interaction. Recently, metallic alloys have been combined to form 
high-temperature ductile/ductile composite materials, raising the natural question of 
whether these metallic composites exhibit the same phenomenological features as do their 
monolithic constituents. Here, this question is addressed in the context of a limited, yet 
definitive (to illustrate creep/plasticity interaction) set of experimental data on the 
model MMC system W/Kanthal. Furthermore, it is demonstrated that a unified 
viscoplastic representation, extended for unidirectional composites and correlated to 
W/Kanthal, can accurately predict the observed longitudinal composite creep/plasticity 
interaction response and strain rate dependency. Finally, the predicted influence of fiber 
orientation on the creep response of W/Kanthal is illustrated. 


NOMENCLATURE 

Stresses 


s ij 

a ij 


Y 

K 


the Cauchy stress tensor 

the deviatoric stress tensor 

the deviatoric internal (or back) stress tensor 

the effective deviatoric stress tensor 

the reference uniaxial creep stress 

tensile threshold stress 
shear threshold stress 


Invariants 

J 2 second invariant of the deviatoric stress (expressed for shear) 

1^ invariant representing transverse shear stress 
I 2 invariant representing longitudinal shear stress 

I invariant representing maximum normal stress in the fiber direction 

o 

F Bingham-Prager threshold function 



Subscripts 


( longitudinal properties 

( transverse properties 

i,j tensor indices, take on values of 1,2 or 3. 

Material Parameters 


/i represents the viscosity of material 
H measure of the hardening 

R measure of thermal recovery 

n,m,/? power law exponents 

E elastic modulus 

v Poisson ratio 

G shear modulus 


Miscellaneous 


dj unit vector denoting local fiber direction 
D.j second order direction tensor; d-dj 
<5.j Kronecker delta function 


€jj inelastic strain rate tensor 

lj ratio of longitudinal to transverse tensile stress 

77 ratio of longitudinal to transverse shear stress 

p smoothing function for the internal stress during cyclic loading 
r time at which load history is constant 

t time 

< > Macauley bracket operator 
Hv[] Heaviside unit function 


INTRODUCTION 

Structural alloys enter the creep range when the service temperature reaches 
approximately 0.3 to 0.4 of their melting temperature, T m . In the terminology of Ashby 
deformation maps [ 1 ], this temperature delimits the boundary of the dislocation creep 
region. At temperatures within the creep range, as at lower temperatures, inelastic 
behavior results from dislocation motion, generation, and interaction but now with time- 
dependent manifestations of diffusion coming into play (e.g., dislocation climb, thermal 
recovery, etc.). 

Aithougn this represents the understanding of the material scientist regarding the 
creep of metals, it only recently has found its way into the mathematical representations 
(constitutive equations} that structural analysts and design engineers use to describe 
material behavior. The recognition that creep and plasticity have a common 
microstructural origin in this temperature range has led to the development of unified 
constitutive models [2,3]. 

The need for unified representations was recognized more than a decade ago in the 
results of phenomenological tests that exhibited strong creep -plasticity interaction [ 2 ]. 
Creep response was observed to be intimately affected by prior history of plasticity and 
vice versa. In terms of the mathematical forms appropriate for describing hereditary 
behavior, these observations prompted the introduction of internal state variables such as 
the internal (or back) stress . In this context the internal stress represents a macroscopic 
measure of the stress field associated with the dislocation microstructure. Incorporating 
(Bailey- Or owan) evolutionary laws that reflect competition between hardening and 
recovery mechanisms, unified constitutive models are found to accurately predict creep, 
rate sensitive plasticity and their mutual interactions. 

When metallic alloys are combined to form a high-temperature ductile/ductile 
composite material, a natural question arises: Do these metallic composites exhibit the 
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same phenomenological features (e.g., creep -plasticity interaction) that led to unified 
constitutive relationships for their monolithic constituents ? The objective of this study is 
to answer this question in the context of a limited, yet definitive (to illustrate 
creep/plasticity interaction) set of experimental data on a model metal matrix composite 
(MMC) system. We shall also demonstrate that a unified viscoplastic representation, 
extended for unidirectional composites, can accurately predict the observed composite 
response. 


EXPERIMENTS 

A general experimental approach to investigating history-dependent behavior is to 
subject identical specimens (over a time period t = 0 to t = r) to loading histories 
that differ in some respect. Then for t > r, each specimen is subjected to the same 
loading conditions, during which time their respective responses are compared. Only if 
the material has no memory of its earlier loading history during t < r will the 
responses be the same for t > r. Well defined tests of this kind are not only useful in 
studying hereditary behavioral features but also in identifying appropriate state variables. 
Here, we focus on a pair of tests of this kind [4] that are designed specifically to 
identify creep/plasticity interaction effects. In tnese tests the different initial histories (t 
< r) of the two identical specimens correspond to different amounts of plasticity on 
load-up. The common loading (t > r) is creep at the same constant stress (cr Q ). 

The test material is a unidirectional W/Kanthal composite of 35 percent fiber 
volume ratio. Plate (coupon) specimens were tested at 600° C under a uniaxial stress 
along the fiber (0‘) direction. This temperature corresponds to approximately 0.2 T ra for 
the tungsten fibers and 0.45 T m for the Kanthal matrix. The load-up and creep 
responses for both of the test specimens are shown in Fig.l; both tests are performed 
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Strain, e, percent 

(b) Stress strain response of 
specimen 2. 

Figure 1. — Experimental observation of creep- plasticity interaction in W/Kanthal at 600 °C with a 
fiber volume ratio of 35%. 
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under load control, with an engineering stress rate of 0.193 MPa/s. In each case, the 
resulting elastic strain rate on load-up is approximately 1 /xc/sec. Specimen 1 is loaded 
directly to the creep stress a Q = 520 MPa; specimen 2 is first overloaded to 610 MPa 
and then unloaded to a 0 = 520 MPa. The subsequent creep responses during the 4.1 
hour hold period, AB for specimen 1 and A'B' for specimen 2, are also shown in 
Fig.l. Evidently, creep varies markedly with load-up history. The initial creep rate of 
specimen 2, whose loading history involves the overload, is at least two orders of 
magnitude less than that of specimen 1; see Fig. lc. A significant permanent strain is 
also shown in each specimen upon unloading. The W/Kanthal composite viewed 
phenomenologically as a material in its own right, evidently exhibits plasticity, creep and 
their mutual interaction, even as tested here along the strong fiber (0*) direction. 


UNIFIED CONSTITUTIVE MODEL 

We now state an anisotropic deformation theory to be used in modeling the 
W/Kanthal behavior of Fig.l. The continuum based theory [5,6] considers the composite 
as pseudo homogeneous and locally transversely isotropic with its own properties that 
can be measured for the composite as a whole. A complete experimental procedure for 
characterizing a particular composite in terms of this theory is specified in [5], An 
isothermal, multiaxial statement of the model is as follows: 

Flow law 


2 M «ti = <F) n r i j 

Evolutionary law 


in which, 


n ii — 


H 


iij - R ij 


K| l 


G = 


K<|> l 


Ii + H — Ij - 1 

II + — jjjl2 + — — — 1 3 

11 4(4o/ 2 -l) J 


r ij = Sij - £(DkiEjk + DjkSki - 2I„Dij) 

— s-CIoCSDjj — tfy) 


^ij — 3-ij “ £(DkiUjk + Djk&ki — 2I 0 Dij) 

— 5CIo(3Dij — iij) 


1 

t 

< 


= < G - G 0> 



Hv[S ijTij J + G o 
0 < £ < 1 
o < c < i 


(1) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 
( 9 ) 


where ui = Y^/Y,p and 77 = K^/Krp represent the ratio of longitudinal to transverse 

tensile and longitudinal to transverse shear threshold stresses, respectively. These 
equations incorporate the following invariants, whose physical meanings are discussed in 

[6J: 
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with 


and 


11 — J2 - I + ^ I3 

12 = I - I3 

1 3 - n 

J2 = ^ijEji 

I = DijEjkSki 

I 0 = DijEji 

Eij — S ij — aij 

Djj = <1 id j 


( 10 ) 

( 11 ) 

(12) 

(13) 

( 14 ) 

(15) 

(16) 
(17) 


where eg is the inelastic strain, Sij is the applied deviatoric stress, a y is the 
deviatoric internal stress and di is a unit vector designating the local fiber direction 

[6]. The invariants 1 1, I2 and I3 in Eqs. (1) to (6) are found from Eqs. (10) to (15) 
with Eji replaced by a§j. 

The material constants are: 


/i, n, H, /?, R, m, K T , G Q , v, and rj 


As stated earlier, they must be determined for a specific composite using the procedures 
outlined in [5], 

A compatible anisotropic elasticity theory is necessary for the subsequent 
correlations and predictions. Here, we use the anisotropic elasticity model developed by 
Arnold [7]. The necessary elastic moduli and Poisson ratios are: 

E l , E,p G l , V L and 

where the subscripts L and T denote longitudinal (along the fiber direction) and 
transverse (normal to the fiber direction), respectively. 


MODEL CHARACTERIZATION - SPECIMEN 1 

Given that the complete set of characterization tests specified in [5] are not 
available for W/Kanthal at 600* C, a limited data set consisting of the tensile and creep 
response of specimen 1 (Figs, la and lc) is used. Characterization is thus accomplished 
by reducing the multiaxial equations of the previous section to the uniaxial (0° ) form 
and then determining optimal values of the material constants to match key features of 
the experimental data, viz., the tensile load-up and the transient creep response of 
specimen 1. The procedure outlined in [5] was followed as closely as possible allowing for 
the very limited data set. The values of the relevant constants so obtained are: 


Inelastic 

Elastic 

[L = 2.75 x 10 13 (3.0/(4w 2 -l)) 
n = 5.3 

E L = 206850. 

H = 3.45 x 10 5 

E t = 158047. 

/? = 0.55 

0 

II 

R = 8.87 x 10' 8 ((4w 2 -l)/3.0) 

u T = 0.25 

m = 1.63 

G l = 68950. 


K| = 4.543(3.0/ (4o/ 2 — 1)) 

G = 0.05 
0 
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These values are consistent with the units of MPa for stress and hr for time. 
As the predictions in the following section are made for the same uniaxial direction (O') 
along which the correlations are based, the material constants u and 77 , designating 
the degree of anisotropy, need not be specified. Figures 2a and 2c show a comparison 
of the experimental (open circles) and predicted responses (dashed lines) of specimen 1 
based on this limited data correlation. The accuracy of the correlation is very good. 



Figure 2. — Comparison of experimental (symbols) and analytical predictions (lines) of creep- 
plasticity interaction in W/Kanthal at 600 °C with a fiber volume ratio of 35%. 


PREDICTED RESPONSE - SPECIMEN 2 

Having used the data of specimen 1 for characterization, we now employ the 
constitutive model with the specified constants to predict the response of specimen 2. 
This is done by numerically integrating the constitutive equations over the prescribed 
loading history. Figures 2b and 2c compare this prediction (solid lines) with the 
experimental response of specimen 2 (open circles). The agreement is remarkably good, 
showing definitively that the unified model [5,6] accurately represents the observed 
hereditary effects and the interaction of plasticity and creep in these experiments. 

The experiments described here were conducted exclusively along the 0°-fiber 
direction of the unidirectional W/Kanthal composite. If inelastic/hereditary effects 
were to be suppressed, it is expected this would occur under testing, as here, along the 
strong (O’ ) fiber direction. Evidently, this is not the case . However, this raises the 
question of how hereditary features might vary as the anisotropic composite is loaded in 
other directions fe.g., normal to the fiber direction). Such predictions for other fiber 
orientations can be made for W/Kanthal at 600' using the results obtained here, 
provided the remaining constants u and 77 are specified [5]. The determination of 
these constants, characterizing the degree of anisotropy, and the required experiments 
analogous to the pair described here but at, say 90’, are left for future research. 
Qualitative predictions of the influence of fiber orientation are made in a subsequent 
section based on assumed values of u and rj. 
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STRAIN RATE DEPENDENT PLASTICITY 


Another feature of elevated temperature response of monolithic metals is strain rate 
dependence. This feature can also be represented by unified models [2,3]. As with 
hereditary effects, we might ask if strain rate dependence, such as observed in monolithic 
alloys, similarly carries over to MMCs. 

To address this question, we include the load-up response of a third W/Kanthal 
specimen tested (under strain control) at 600° C and at a strain rate of 100/x£/sec, (two 
orders of magnitude greater than that for specimens 1 and 2, described earlier). Figure 3 
compares the tensile response of this specimen with that of specimen 2. Rate 
dependence is evident , showing approximately a 15-20% increase in strength for the 
100-fold increase in strain rate. This is approximately the same degree of rate 
dependence exhibited by many monolithic alloys at temperatures in this range. 



Figure 3.— Comparison of experimental (symbols) and analytical pre- 
dictions (solid lines) of stress-strain response at different strain rates 
of 1 ^e/sec and 1 00 p-t/sec. 


INFLUENCE OF FIBER ORIENTATION 

Finally, the predicted influence of fiber orientation on the creep response of 
W/Kanthal is illustrated in Fig. 4 for fiber angles (relative to the loading direction) of 
0, 30, 60 and 90°. These calculations are made with assumed values of u = 5.0 and 77 
= 2.5 designating the degree of anisotropy and so are to be considered only qualitative. 
The creep stress in Fig. 4 is 200 MPa and the equivalent elastic strain rate is 1/u/sec. 

As expected, the inelastic response significantly increases, during load-up as the 
fiber angle is increased, see Fig. 4a, such that in the transverse case (90° ) a total strain 
equal to approximately 0.4% is achieved at 200 MPa instead of 520 MPa, as in the 
longitudinal (0° ) case shown in Fig. 2a. This difference in maximum achieved stress 
(by a factor of 2.5) indicates a significant decrease in transverse load carrying capacity 
relative to that in the longitudinal direction, as expected. In comparing the 60° case to 
that of the 90° case (see, Fig. 4b) it is interesting to note that the amount of creep 
strain incurred in the 60’ case exceeds that incurred in the 90° case, whereas the 0° 
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and 30* cases are ordered as expected. Past experience has shown, however, that this 
change in ordering is dependent upon the relative magnitude of the material parameters 
chosen. Clearly, if u and f? (the degree of anisotropy] are found experimentally to differ 
greatly from the above assumed values, so too will the predicted amount of inelasticity. 



90 

60 



Time, hr 

Figure 4.— Analytical predictions of the creep response of W/Kanthal at varying fiber angles, e : 
0, 30, 60, 90°, with a fiber volume ratio of 35%. 


CONCLUSIONS 

Generally, we conclude from this study that at 600“ C W/Kanthal, viewed here as 
a model ductile/ductile MMC, shows qualitatively the same features of time-dependent, 
hereditary behavior observed in monolithic alloys at comparable temperatures. That is to 
say, WyKanthal exhibits rate sensitive plasticity, creep and their intimate interaction. 
These effects are evident even as tested here uniaxially (along the direction of maximal 
constraint, i.e., 0“ ) in the fiber direction. 

We further conclude that the unified anisotropic viscoplasticity theory [6] used to 
model the W/Kanthal response is capable of representing the key features of the 
inelastic behavior, for example creep-plasticity interaction, strain-rate sensitivity, and the 
effect of fiber orientation. Here, we have demonstrated, with a limited set of data, the 
accuracy of this model in the longitudinal direction and leave for future work the 
validation of the off-axis response (i.e., effect of fiber orientation). 
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